Ground state of spin-1 Bose-Einstein condensates with spin-orbit coupling in a 
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We systematically investigate the weakly trapped spin-1 Bose-Einstein condensates with spin- 
orbit coupling in an external Zeeman field. We find that the mean-field ground state favors either 
a magnetized standing wave phase or plane wave phase when the strength of Zeeman field is below 
a critical value related to the strength of spin-orbit coupling. Zeeman field can induce the phase 
transition between standing wave and plane wave phases, and we determine the phase boundary 
analytically and numerically. The magnetization of these two phases responds to the external 
magnetic field in a very unique manner, the linear Zeeman effect magnetizes the standing wave 
phase along the direction of the magnetic field, but the quadratic one demagnetizes the plane wave 
phase. When the strength of Zeeman field surpasses the critical value, the system is completely 
polarized to a ferromagnetic state or polar state with zero momentum. 

PACS numbers: 05.30.Jp, 03.75.Mn, 67.85. Fg, 67.85. Jk 
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I. INTRODUCTION 

Spin-orbit coupling, which is generally referred to the 
coupling of a particle's spin with its degree of motion 
in quantum physics, is responsible for many fundamen- 
tal physical phenomena in quantum systems, such as the 
spin Hall effects and topological insulators Recently, 
with the pioneering experimental realization of an arti- 
ficial abelian or non-abelian gauge potential in neutral 
atoms @, the effective spin-orbit coupling has been 
created in spinor Bose-Einstein condensates (BEC) by 
dressing two atomic spin states with a pair of lasers [J- 
0, |1| , and has attracted a great deal of attention in the 
condensed matter community. Especially, the spin-orbit 
coupling effects can give rise to many intriguing exotic 
ground states in an interacting ultracold spinor Bose gas, 
such as the plane wave (PW), standing wave (SW), tri- 
angular lattice and square lattice phases |9l— f29jj . More- 
over, in pseudo spin-1/2 BEC, the combined effects of 
Rashba spin-orbit coupling and rotating trap can also 
produce some unusual topological patterns including the 
skyrmion and giant vortex [30|-|32|. 

Generally, without spin-orbit coupling, a spinor BEC 
exhibits a variety of magnetic phenomena. Taking the 
spin-1 BEC as an example, the mean- field ground state 
can be a ferromagnetic, antiferromagnetic or polar states 
;33J . Which type of these phases is favored by the ground 
state depends on the spin-dependent atomic interaction. 
Interestingly, an external Zeeman field can transform 
these phases into each other via tuning their magneti- 
zation, and results in a rich ground state phase diagram 
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[3J, |35| . For example, a broken-axisymmetry state can 
emerge as an intermediate phase in the transition be- 
tween the ferromagnetic and polar phases. Such phase 
transitions reflect that the magnetic system responds to 
the external magnetic field in a very unique manner. So 
the application of an external magnetic field provides a 
powerful way to experimentally manipulate the magnetic 
behaviors of a spinor BEC [34{ . 

However, when the spin-orbit coupling is taken into 
account, the situation changes. The ground state of 
spin-1 BEC now could be a PW phase occupying a sin- 
gle momentum state or a SW phase formed by the co- 
herent superposition of two plane waves with opposite 
momentum 10], and these two phases possess different 
magnetic properties. So it is interesting and impor- 
tant to study the effects of an external Zeeman field on 
these phases, which is the main context of this work. 
To make this point more transparent, in this paper, we 
consider a weakly trapped two-dimensional spin-1 BEC 
with Rashba spin-orbit coupling in a Zeeman field. We 
find that, if the Zeeman field is not too strong, the sys- 
tem is in a magnetized SW phase or PW phase, and 
the phase boundary immerses into SW phase due to the 
competition between Zeeman effects and spin-dependent 
atomic interaction. The Zeeman field modulates the 
magnetization of these two phases in different manners. 
For example, the linear Zeeman field magnetizes the 
SW phase, but the quadratic one demagnetizes the PW 
phase. While the strength of the Zeeman field surpasses 
a critical value determined by the strength of spin-orbit 
coupling, the system is completely polarized to either a 
ferromagnetic state or a polar state with zero momen- 
tum. 

This paper is organized as follows. We first give the 
model in Sec. II, and discuss the single-particle ground 
state in Sec. III. By using the variational approximation 
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method and the numerical simulation, we give the phase 
diagram of mean-field ground state and study the effects 
of Zeeman field on the magnetization of this system in 
Sec. IV. Finally, we conclude our results in Sec. V. 



II. THE MODEL 

We consider a quasi- two-dimensional spin- 1 BEC with 
Rashba type spin-orbit coupling in the Zeeman field, 
where the particle's spin couples its degree of motion in 
xy plane. In the mean-field approximation, the Gross- 
Pitaevskii energy functional of such system is of the form 
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+ j-f [ipi {-id x - d y )tjj +Vo {-idx - d y )ip-i + h.c] 

+ 7^+7 IK " «-i) 2 +2|^^o + ^o>-il 2 ]}, (i) 

where m = 1,0, —I and the density distribution of the 
component m reads n m — |^ m | with the condensa- 
tion wave function ip m , thus the total atomic density is 
n = X)m"™- For the spin-orbit coupling term, we con- 
sider the symmetric Rashba case k x — n y = k, where 
k is the strength of spin-orbit coupling. For the inter- 
action terms, c$ and ci are the interaction parameters 
which depend on two-body s-wave scattering lengths do 
and a 2 for total spin 0,2: Co = ATrh 2 N(ao + 2a2)/3M 
and C2 = Anh 2 (a2 — ao)/3M with the atomic mass M 
and total atomic number N. To avoid the collapse of 
BEC under attractive interaction, cq is assumed to be 
non-negative throughout this paper. Moreover, p and q 
represent the strength of linear and quadratic Zeeman ef- 
fects, respectively. For the sake of simplicity, the external 
magnetic field is assumed to be applied in the z-direction, 
and the values of p and q are taken to be non-negative 
constants. 



III. SINGLE-PARTICLE GROUND STATE 

To gain some intuitions, it is instructive to start our in- 
vestigations with the noninteracting Hamiltonian in the 
homogeneous case. The single-particle Hamiltonian pre- 
serves the symmetries of simultaneous spin and space ro- 
tations around the direction of applied magnetic field. 
The group describing such symmetries is 50(2) which 
has discrete subgroups C nz [III- The single-particle 
ground state in absence of the external magnetic field 
is not unique so that it is infinitely degenerate along a 

\Z2k/H in 

where k x and 
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circular ring with radius k± = ^Jk 2 

two-dimensional momentum space [13, \22 
k y are the components of momentum along x and y di- 
rections, respectively. In presence of the Zeeman field, 
we first consider two simple cases in which the linear and 




FIG. 1: (Color online) (a) The single-particle energy spectrum 
in linear Zeeman field with p = 2. (b) The single-particle 
energy spectrum in quadratic Zeeman field with 5 = 1. Noting 
that the figures in the bottom of (a) and (b) represent the 
projections of E- on the k x -k y plane, and the white circular 
rings are the momentum of infinitely degenerate ground states 
in E-. Furthermore, to clarify the energy gap between these 
energy branches, figures (c) and (d) have plotted the energy 
spectrums at k x = 0. In all the figures, the other parameters 
are the same that M = H = 1 and re = 2. 



quadratic Zeeman effects are considered separately. The 
combinative effects of them will be considered later. 

Due to the typical Rashba spin-orbit coupling, the 
single-particle energy spectrums in momentum space can 
be obtained by virtue of the exact diagonalization scheme 
as follows 



L: E = E k , E ± =E k ± X p 
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where L and Q represent the energy spectrums in lin- 
ear and quadratic Zeeman fields, respectively. E^ — 
h 2 k 2 _/2M is the kinetic energy of a free particle. The 
energy spectrums in Eq. ^ have been shown in Fig. 
[2a) and [Ub) , we can observe that the energy branch 
_E_ has lower energy no matter in the linear or quadratic 
Zeeman field. By minimizing the energy branch E_ with 
respect to fcj^, we can obtain the momentum of single- 
particle ground state which are kj^ = ij 20. _ and 

k± = *\/32k 4 — 2q 2 M 2 / (4hn) , respectively. The corre- 
sponding single-particle eigenstates are given by 
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FIG. 2: (Color online) (a) The momentum of single-particle 
ground state changes with linear and quadratic Zeeman ef- 
fects, (b) The coefficients a m of single-particle eigenstates in 
Eq. (J3]) change with the linear and quadratic Zeeman effects. 
In both (a) and (b), the parameters are k = 1 and h = M = 1. 



where 9 = arctan (k y /k x ), and the coefficients a m satis- 
fying Em "m = 1 ar e 

2 K 2 + Mp 2k 2 -Mp r- 

1 ai= — 4^2 — ' a - 1= — 7^2 — I «o=V 2Q; i Q! -iA 4a ) 

^4 K 2 - Mq V2^/4 K 2 + Mg 

Q : ai=a_i= ,a = -; . (4b) 
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From the results above, on the one hand, the Zeeman 
effects shift the momentum of ground state compared 
to the case in absence of external magnetic field. All the 
single-particle states with same momentum k± but differ- 
ent azimuthal angle in two-dimensional momentum space 
are degenerate ground states, namely, the single-particle 
ground states are infinitely degenerate along a circular 
ring with the radius k±_ shown in Fig. HJa) and[T^b). 
Moreover, since the momentum of ground state should 
be real, we have noted that there exists two additional 
restraint conditions p < 2k 2 /M and q < 4k 2 /M. In 
particular, if the strength of linear or quadratic Zeeman 
effect takes the maximal critical value, the circular ring of 
momentum in infinitely degenerate ground states shrinks 
to the zero momentum point (the center of ring shown 
in Fig. (Ha) andQJb)) so that the system is polarized to 
a non-degenerate ground state. On the other hand, the 
presence of Zeeman field can open a energy gap between 
these energy branches as shown in Fig. [IJc) andQJd). 
In linear Zeeman field, the gap between adjacent energy 
branches are A + o = E+ — Eq = p and Ao- = Eq — E_ = 
K 2 /M+Mp 2 /4n 2 . While the quadratic Zeeman field only 
opens the gap A+_ = A _ = q/2 + k 2 /M + Mq 2 /l&n 2 
between the energy branches E+ (or Eq) and E-, it is 
gapless between the energy branches E + and Eq. These 
facts demonstrate that the linear and quadratic Zeeman 
effects play different roles on affecting the structure of 
single-particle energy spectrum. 

Let us now consider the combination effect of the lin- 
ear and quadratic Zeeman terms. For simplification, we 
just give the lower energy branch in single-particle energy 



spectrums as follows 

E- = E k + ^ +uj\fVA-<d + cj 2 \f-VA-e, (5) 

where oj = - exp (iir/3), A = 6 2 - T 3 with 9 = q(q 2 + 
18K 2 E w /m ~ 9p 2 )/27 and T = {3p 2 + q 2 + I2n 2 E^/m) /9. 
We have to point out that p and q must satisfy A < 
for the real energy eignvalue E-. Since the expression 
of energy spectrum in Eq. ([5]) is very complex, it is im- 
possible to analytically solve the momentum of single- 
particle ground state and derive the coefficient a m of 
the corresponding eigenstates in Eq. ([3]), so we plot this 
quantities in Fig. [5] by seeking the minimal value of en- 
ergy Eq. (|5|) numerically. It can be clearly seen that 
the quadratic Zeeman effect dominates in region I, but 
the linear Zeeman effect dominates in region III. In both 
the regions I and III, the momentum of single-particle 
ground state is zero, the corresponding coefficients a m 
in single-particle eigenstate are (0,1, 0) T and (1,0, 0) T , 
respectively. These results imply that the system is po- 
larized by the magnetic field in both the regions I and III. 
While in region II, the coefficient satisfies < a m < 1, 
and the momentum of ground state is nonzero so that 
the single-particle ground states are infinitely degenerate 
along a circular ring with radius k± in momentum space. 
Furthermore, the boundary between the regions II and 
III can be analytically fixed that isq/3— p — ( = with 

C = u V7a -o + lj 2 V-Va - e. 



IV. MEAN-FIELD GROUND STATE IN 
EXTERNAL ZEEMAN FIELD 

For the spin-1 BEC with Rashba spin-orbit coupling, 
the interaction effects in absence of Zeeman field have 
been investigated extensively (§-(2§|. In homogeneous 
case, the mean-field ground state favors either the PW 
for C2 < or SW for ci > 0, the former is ferromagnetic 
state, but the latter is polar state [Toj ] . These phases 
spontaneously break the rotational symmetry around the 
z direction despite the fact that the Hamiltonian is ax- 
isymmetric. In particular, these two different phases are 
degenerate for C2 = 0. Now the point is that if a Zeeman 
field along z-direction is turned on, we are inquisitive 
about whether an external magnetic field can essentially 
change the phase diagram of mean-field ground state for 
spin-orbit coupled spin-1 BEC? How the external mag- 
netic field affects the magnetization of ground state? The 
investigation of the unique features of these problems is 
the primary purpose of this section. 

In two dimensional homogeneous system, the inter- 
acting Hamiltonian preserves the symmetries of single- 
particle Hamiltonian, and the atomic interactions can 
couple different single-particle states so that the ground 
state preserving symmetry C nz can be approximately 
described by the linear superposition of single-particle 
eigenstates on the degenerate momentum ring, \& = 
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FIG. 3: (Color online) (a) The phase diagram of mean-field 
ground state in linear Zeeman field, (b) The phase diagram of 
mean-field in quadratic Zeeman field. The red triangles and 
green squares respectively denote the SW and PW phases 
given by numerical results, and the blue line is the phase 
boundary £ = between these two phases. The blue circles 
represent the zero momentum ferromagnetic state (FS) in (a) 
and the zero momentum polar state (PS) in (b), respectively. 
The other parameters in both the figures are Co = 100, M = 1 
and k = 6. 

E„A*k with Ai satisfyin g Y \A n \ 2 = 1, where n 
is an non- negative integer [ljj, |26| . By substituting the 
ansatz into Eq. ((TJ) to compute the energy, we find that 
these high-symmetry states with n > 3 are energetically 
unfavored by ground state which can be verified by our 
numerical results. Thus this fact motivates us to choose 
the following ansatz for clarifying the phase diagram of 
mean-field ground state and developing a simple physical 
understanding in below discussions 

* = -4* k + £5<_ k , (6) 

where \l/k and \I/_k are two counter-propagating plane 
waves in Eq. (O, and the superposition coefficients A 
and B are real constants satisfying the normalization con- 
dition A 2 + B 2 = 1 under £ m / |V,„| 2 dr = 1. 

By inserting the ansatz into the mean-field energy 
functional Eq. ([T]), we can obtain the energy expression 
as a function of the parameter A 

E[A] =£(A i - A 2 ) + const, (7) 

where £ = [2al(al + 6ai«_i + a 2 _ x ) - (a? - a!_ 1 ) 2 ]c 2 - 
(a 2 — 6Vq + oPl-^cq. Noting that some constants unre- 
lated to A in Eq. ([7]) are omitted. We can immediately 
determine the phase diagram of mean-field ground state 
by minimizing the energy with respect to A: (i) The PW 
phase with A = and B = 1 if £ < 0; (ii) The SW phase 
with A = B = l/v2 if £ > 0. In particular, these two 
phases are degenerate when £ = 0. These results demon- 
strate that the external magnetic field shifts the phase 
boundary from c 2 = to £(c 2 ,p, q) = due to the com- 
petition between the Zeeman effects and spin-dependent 
interaction. 

To further explore the physics behind this, in what fol- 
lows, we perform the detailed analysis on the phase dia- 



gram of ground state, and investigate the effects of exter- 
nal magnetic field on the magnetic behaviors of this sys- 
tem. We first consider the linear and quadratic Zeeman 
effects, respectively. The combination effects of them will 
be considered afterwards. 



A. Phase diagram and magnetic behaviors of 
ground state in linear or quadratic Zeeman field 

In the linear or quadratic Zeeman field, the phase 
boundary are £ = (2 - M 2 p 2 /n 4 ) c 2 + (c 2 - c ) or 

£ = 2c 2 — -j^Si (2c2 + c ), respectively. In parameter 
spaces of (c2,p) and (c2,q), the corresponding phase di- 
agrams of mean-field ground state have been shown in 
Fig. |31 We can observe that the SW phase only ex- 
ists in the case of c 2 > 0, and the region of which is 
less than the region occupied by the PW in the c%-p 
or C2~q plane. More importantly, the external magnetic 
field can induce the phase transition from the SW to PW 
when the strength of magnetic field surpasses the criti- 
cal value p a = 2k 2 ^ ( v / 2q[+ 2c 2 c - 2c 2 )/[M 2 (c - c 2 )] 
on the phase boundary for the linear Zeeman effect (or 
qo = 4k 2 -\/2c 2 /[(co + 2c 2 )M 2 ] for the quadratic one), 
where we require cq 3> c 2 in accordance with the real 
experiments. In particular, if the strength of linear 
(or quadratic) Zeeman effect attain the critical values 
p m = 2k 2 I'M (or q m = 4k 2 /M), the system can be po- 
larized by the magnetic field so that all the atoms will 
occupy the state ipi (or tpo) with zero momentum, in 
which the ground state is ferromagnetic state (or polar 
state). 

Based on the coupled Gross-Pitaevskii equations, we 
have checked all the analytical predictions by using the 
imaginary time evolution method to numerically look for 
the ground state solutions which can minimize the mean- 
field energy. We would like to point out that a very weak 
isotropic harmonic potential is included in our numerical 
simulation, it does not change the results in the homo- 
geneous case. As shown in Fig. [21 these analytical pre- 
dictions given by the ansatz in Eq. ([6]) agree with the 
numerical results very well, which demonstrates that the 
excellent ansatz in Eq. © captures all the fundamental 
physics. 

In hindsight, the shifts of phase boundary can 
be understood by analyzing the magnetization of 
ground state in response to the external magnetic 
field. Generally, the magnetization for spin-1 BEC 

is defined as M = J ' M 2 + M 2 + M 2 Z with = 

J2mnl 1 Pm( cr n)mn4>ndxdy being the component of mag- 
netization along the directions [i = x, y, z, where the sub- 
scripts m, n — 1, 0, — 1, and represents the 3x3 spin 
matrix for spin-1 (3o| . A priori, we can use the ansatz in 
Eq. ^ to analytically compute the magnitude of mag- 
netization for the cases of c 2 < and c 2 > 0, respec- 
tively. The results are listed in the table It can be 
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Linear Zeeman field 


Quadratic Zeeman field 


C 2 > 


A/fP<Po — M P 
JVl SW — 2^ 




KAP>VO _ 1 

J vl pw — 1 


JVl PW - iK i 


C 2 < 


Mpw = 1 


M.pw = -* 



TABLE I: The expressions of magnetization in linear and 
quadratic Zeeman fields, respectively. For C2 > 0, the ground 
state during p < po or q < go is the SW phase, but the PW 
phase for p > po or q > qo. For ci < 0, the ground state is 
always the PW for an arbitrary p or q. 



seen that the SW and PW phases respond to the exter- 
nal magnetic field in very different manners. On the one 
hand, in linear Zeeman field, the magnitude of magne- 
tization of PW phase for C2 < always keeps invariant 
regardless of the external magnetic field. While in the 
regions p < po and ci > 0, the system is initially in 
SW phase, the magnitude of magnetization increases lin- 
early from zero until p > po, in which the system has 
transformed into PW and the magnitude of magnetiza- 
tion undergoes a discontinuous sudden transition from 
M = Mp /(2K 2 ) to M. = 1. On the other hand, an 
opposite story occurs in quadratic Zeeman field. With 
increasing the strength of magnetic field, the magnitude 
of magnetization of PW phase for C2 < decays quadrat- 
ically from A4 = 1 to A4 = 0, which implies that the 
quadratic Zeeman effect demagnetizes the PW. However, 
in SW phase for q < qq and ci > 0, the quadratic Zee- 
man effect can not change its magnetization with Ai = 0. 
When the strength of magnetic field surpasses the critical 
value go, the system s uddenly possesses magnetization 
M = ^/16k 4 — M 2 Qq/(4:K 2 ), which demonstrates that 
the phase transition from SW to PW occurs. Again, 
we have verified these predictions numerically and the 
results are summarized in Fig. |4l 

To highlight the effects of Zeeman field on the mag- 
netization of ground state, further, we map the magne- 
tization vectors of PW and SW phases onto a sphere 
with radius M. — 1 as shown in Fig. [SJ Before perform- 
ing the analysis in detail, we introduce the polar angle 
d = arctan ^ and azimuth angle tp = arctan of 
the magnetization vector in spherical coordinate frame, 
where A4+ = A4 X + iM. v . In presence of the external 
Zeeman field, the polar and azimuth angles of the mag- 
netization vector for these two phases are 



PW : = arctan ■ 



if : 



(8a) 



a.\ — a~i 

SW:-d = ip = undefined for p = or < q < <?o,(8b) 
i9 = 0, (p = undefined for < p < po, 

7T 

V = -,<p = for q > q Q , 

where 9 = arctan (k y /k x ) is determined by the momen- 
tum of ground state. Motivated by the numerical results, 
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FIG. 4: (Color online) (a) The magnetization and atomic 
distribution of the ground state change with the strength of 
linear Zeeman field in the cases of ci = 10 and C2 = —10. 
(b) The magnetization and atomic distribution of the ground 
state change with the strength of quadratic Zeeman field in 
the cases of C2 = 10 and C2 = —10. In both (a) and (b), the 
solid red and blue lines are respectively the corresponding 
magnetization of C2 > and C2 < given by the analytical 
results shown in table (P), the green lines represent the atomic 
numbers N m . The red stars and blue circles given by our 
numerical results are the magnetization for C2 > and C2 < 0, 
respectively. Moreover, the other parameters in both (a) and 
(b) are the same that k = 6, h — M — 1 and Co = 100. 



we can fix 9 = ? in the following discussions. Noting that 
since the SW phase has no magnetization when p = 
or < q < qo, so the polar and azimuth angles of the 
magnetization vector for SW phase are undefined in the 
center of sphere. Furthermore, as we shall see below, the 
magnetization of SW phase is along the z-direction when 
< p < po in linear Zeeman field, so the azimuth angle 
of SW phase is undefined too. Fig. [S] shows clearly that 
how the external magnetic field changes the magnetiza- 
tion of two phases. For the case of C2 < 0, the ground 
state is PW, the linear Zeeman effect rotates the mag- 
netization vector from transverse direction (xy plane) to 
longitudinal direction (z direction) along the meridian 
with azimuth angle if = 7r/4, but the length of vector is 
invariant shown in Fig. [SJa). On the contrary, along the 
azimuth angle <p = 7r/4 direction, the quadratic Zeeman 
effect demagnetizes the PW in xy plane by reducing the 
length of magnetization vector from one to zero as shown 
in Fig. [Sfc). For the case of C2 > summarized in Fig. 
[5Jb) andlDJd), the system would be initially in SW phase 
if p < po or q < qo . In this case, the length of magnetiza- 
tion vector in linear Zeeman field increases linearly along 
z direction from zero, but it keeps invariant in quadratic 
one with M. — 0. While p (or q) surpasses the critical 
value po (or go), the SW changes into the PW so that 
the magnetization vector responds to the Zeeman fields 
in the same manner as the one in PW. 
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(a): c 2 < (b): c 2 > 




(c): c 2 < (d): c 2 > 



FIG. 5: (Color online) The magnetization vector responds 
to the linear (L) or quadratic (Q) Zeeman field. In (a)-(d), 
the red and blue solid circles represent the end of vectors for 
C2 > and C2 < 0, respectively. Noting that the direction of 
the magnetization vector points to the red or blue circle from 
the center of sphere. Furthermore, the red and blue solid 
lines are the change trajectories of the end of magnetization 
vectors for SW and PW phases, respectively. In linear Zeeman 
field, the end of vector for C2 > undergoes a transition 
from the z direction to the blue line when p > po- While in 
quadratic Zeeman field, the end of the magnetization vector 
for C2 > undergoes the transition from the center of sphere 
to the azimuth angle direction when q > qo- 



B. Phase diagram and magnetic behaviors of 
ground state in the combination of linear and 
quadratic Zeeman fields 

So far, we have focused on the simple case that the 
linear and quadratic Zeeman effects are dealt with sep- 
arately. Experimentally, the linear and quadratic Zee- 
man effects always coincide with the presence of exter- 
nal magnetic field. The former can be effectively tuned 
by changing the total spin of the system, and the latter 
can be usually tuned by using a linearly polarized mi- 
crowave field due to the alternating current Stark shift 
[3rj H3| • As shown in Refs. [H, HH] , under a certain range 
of linear and quadratic Zeeman effects, the ground state 
phase diagram of spin-1 BEC without spin-orbit coupling 
becomes much richer due to the competition between 
the linear and quadratic Zeeman effects and the spin- 
dependent interaction. The typical phase is the presence 
of a broken-axisymmetry state. Therefore, in this sub- 
section we investigate the phase diagram of spin-1 BEC 
with spin-orbit coupling in the combination of linear and 
quadratic Zeeman effects. 

We now have a similar situation as the single-particle 
problem in this combination effect. Since the expression 
of £ is very complex, it is impossible to analytically fix the 




P 



FIG. 6: (Color online) The phase diagram and magnetization 
distribution in the combination effects of linear and quadratic 
Zeeman fields. In the phase diagram (a), the red triangles, 
blue circles, green squares and purple stars respectively de- 
note the SW, zero momentum ferromagnetic state (FS), PW 
and zero momentum polar state (PS) given by numerical re- 
sults. The dark line is the phase boundary £ = between 
these phases, (b) the magnetization distribution of the sys- 
tem versus p and q. In both the figures, the parameters are 
Co = 100, C2 = 10 and k = 6. 



phase boundary between the SW and PW phases. So we 
plot the phase diagram and magnetization distribution 
numerically in Fig. [5] based on the ansatz Eq. It 
should be pointed out that for the case of c 2 < 0, the 
ground state is always the PW in the whole parameter 
space (p, q) due to £ < 0. Thus we just consider the 
case of c 2 > in following discussions for the sake of 
simplicity. From the phase diagram in Fig. [BJa), we can 
see that the ground state initially favors the SW phase 
when the strengths of linear and quadratic Zeeman effects 
are relative weak. However, the PW phase dominates if 
the quadratic Zeeman effect goes far beyond the linear 
one. Importantly, our numerical results show that all the 
atoms are polarized to the state ipi for q < 14(p — 39)/ll 
with p € [39, p m ] or the state ipo for p = and q = 
q m , respectively. The corresponding order parameters 
of ground state in these regions are = (1,0, 0) T and 
"Jo = (0, 1, 0) T . In these states, the momentum of ground 
state is zero. All the results above can also be equally 
reflected by analyzing the magnetization of the system as 
shown in Fig. EJb), which shows that the two states 
and are ferromagnetic state with M. = 1 and polar 
state with M. = 0, respectively. 



V. CONCLUSIONS 

Within the framework of the mean-field theory, we 
have systematically investigated the weakly trapped spin- 
1 Bose-Einstein condensates with spin-orbit coupling in 
an external Zeeman field, and clarified the effects of an 
external magnetic field on the system. In the single- 
particle case, the linear and quadratic Zeeman effects 
play the different roles on affecting the structure of single- 
particle energy spectrum. When the atomic interactions 
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are taken into account, we give the phase diagram of 
mean-field ground state which is mainly comprised of 
the magnetized standing wave phase and plane wave 
phase. To develop a physical understanding, we use 
the linear superposition of two single-particle eigensates 
with counter-propagating wave vectors as the variational 
ansatz to analytically determine the phase boundary be- 
tween these two phases, which agrees with our numeri- 
cal results very well. The Zeeman field can induce the 
phase transition between the standing wave and plane 
wave phases, such phase transition can be further un- 
derstood by analyzing the response of magnetization of 
these two phases to the external magnetic field. In par- 
ticular, when the strength of Zeeman field surpasses a 



critical value which is related to the strength of spin- 
orbit coupling, the system is completely polarized to a 
ferromagnetic state or a polar state with zero momen- 
tum. These investigations not only help us to deepen the 
understanding of the physics behind the interaction be- 
tween a matter field and gauge field, but also provide an 
effective way to manipulate the spin-orbit coupled spinor 
BEC experimentally by using an external magnetic field. 
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